Charged-current neutrino-^'^^Pfe reactions 
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We present theoretical results on the non flux-averaged Ph{ve,e~)'^°^ Bi and Pb(u^, Bi 
reaction cross sections, obtained within the charge-exchange Random-Phase- Approximation. A 
detailed knowledge of these cross sections is important in different contexts. In particular, it is 
necessary to assess the possibility of using lead as a detector in future experiments on supernova 
neutrinos, such as OMNIS and LAND, and eventually detect neutrino oscillation signals by exploiting 
the spectroscopic properties of ^"^Bi. We discuss the present status on the theoretical predictions 
of the reaction cross sections. 



The study of reactions induced by neutrinos on nu- 
clei is at present an active field of research. A detailed 
knowledge of the reaction cross sections is interesting for 
different domains, going from high energy physics to as- 
trophysics [Q. For example, they are necessary in the in- 
terpretation of current experiments on neutrinos as well 
as in the evaluation of possible new detectors for future 
experiments. The importance of neutrino-nuclei reac- 
tions in astrophysical processes, such as the r-process 
nucleosynthesis, is also being attentively studied [§,1). 
In particular, v — Pb reactions have attracted much in- 
terest recently. Lead has been used as a shielding ma- 
terial in the recent experiments on neutrino oscillations 
performed by the LSND collaboration so that esti- 
mates of the v — Pb reaction cross sections are necessary 
for the evaluation of backgrounds in these experiments; 
also projects on lead-based detectors such as OM- 
NIS @g and LAND §, are being studied for the pur- 
pose of detecting supernova neutrinos. These detectors 
might provide information on neutrino properties, such 
as oscillations in matter or the mass by mea- 
suring the time delay and/or spreading in the neutrino 
signal [HI] as well as help in testing supernova models. 
From the practical point of view, lead-based detectors 
seem to present several of the characteristics required to 
be supernova observatories, namely high sensitivity to 
neutrinos of all flavors, simplicity, reliability with inex- 
pensive materials |9|. Large cross sections for neutri- 
nos in the supernova energy range are also an impor- 
tant condition since they determine the possible rates 
and therefore the maximum observable distance. Actu- 
ally, i/-nucleus reaction cross sections increase strongly 
with the charge of the nucleus. For example, if the neu- 
trinos come from the Decay- At-Rest (DAR) of the 
cross sections of the flux-averaged charged-current (CC) 
reaction +z Xn —*z+i -'^at-i + e~ goes from about 
14. 10-'*2 cm2 for ^^C to 2.56 10 

and is estimated to be 3.62 10 ^ ^ 

sides these practical features which are essential in the 
choice of the nucleus to use to detect neutrinos, another 



important feature is the spectroscopic properties which 
may suggest attractive signals of supernova neutrino os- 
cillations. In ||l^ , for example, it has been shown that the 
measurement of events where two neutrons are emitted 
by '^°^Bi excited in the reaction i/e +^°^ Pb Bi + e" 
is both flavor-specific and very sensitive to the mean en- 
ergy of the Ve- In case when I'p, i^r ^ i^e oscillations take 
place, the hotter v^, would increase the number of two 
neutron events by a factor of fort y [|lO| . Another pos- 
sible signal has been proposed in |l5t], that is that the 
energy distribution of the neutrons emitted in the same 
CC reaction should have a peak at low energy more or 
less pronounced according to whether the oscillations oc- 
cur or not. This peak would come from the excitation of 
a peak at around 8 MeV in the Gamow- Teller strength 
distribution. (One should however note that this peak 
has never been observed experimentally). Both the esti- 
mate of the CC V — Pb reaction cross section in ||l^ and 
the microscopic calculations of ||l5| show that a possible 
oscillation signal relies strongly on the knowledge of the 
spectral properties of ^^^Bi. In fact, the CC reaction 
cross section induced by scales almost as the square of 
the electron energy and is particularly sensitive to the de- 
tailed structure of the excitation spectrum as was already 
pointed out for the case of ^^C ||l6|. It is then important 
either to get the cross sections directly from the experi- 
ment or/and to obtain different theoretical estimates in 
order to know the theoretical uncertainties and how they 
affect the reaction cross sections. This is crucial when the 
impinging neutrino energy increases because not only the 
allowed Gamow- Teller (GT) and Isobaric Analogue State 
(IAS) contribute significantly to these cross sections but 
also forbidden transitions, of first-, second-, third-order 
(which are not very well known experimentally). 

In this paper, we present new theoretical results for 
the CC Ve +^°^ Pb — >208 _|_ g- reaction cross sec- 
tion. Our calculations, as opposed to |ll5|] , are performed 
in a self-consistent charge-exchange Random-Phase- 
Approximation (RPA) with effective Skyrme forces. Con- 
trary to all the previously published calculations, we 
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present non flux-averaged cross sections, obtained for 
both low-energy Ve and high-energy v^. These reaction 
cross sections given as a function of neutrino energy span 
a large energy range. They can be used to convolute with 
different neutrino fluxes in various contexts, for example 
for future experiments with astronomical neutrinos which 
are at present under study, for the very recent terrestrial 
experiments such as the LSND ones ||l5[ to estimate the 
background, or in the r-process nucleosynthesis. 

We will emphasize the importance of the contribution 
of forbidden transitions and how it evolves as a function 
of neutrino energy. This is often not taken into account in 
many present r-process nucleosynthesis calculations and 
so the neutrino-nuclei cross sections are underestimated 
(in jl^ only the importance of flrst forbidden transitions 
in neutron-rich nuclei was emphasized). 

We will compare our results with presently available 
calculations |l0|jr^ ]. With this aim, we will present two 
different flux-averaged cross sections, where the neutrino 
fluxes are given by either the DAR of fX^ and Decay-In- 
Flight (DIP) of 7r+; or by a Fermi-Dirac spectrum for a 
supernovae explosion. Finally, we will discuss our results 
in relation to the suggested possible oscillation signals 
that would use the spectroscopic properties of "^^^Bi. 

The general expression for the differential cross section 
as a function of the incident neutrino energy E,j for the 
reaction n +'^°^ Pb I +2"^ Bi {I = e, n) is @ 

a{E,) = —cos^OcY^PiEi d{cose)Mp, (1) 

TT ^ J-1 

where G cos 9c is the weak coupling constant, 9 is the 
angle between the directions of the incident neutrino and 
the outgoing lepton, Ei ^ Ey ~ Efi (pi) is the outgoing 
lepton energy (momentum), Efi being the energy trans- 
ferred to the nucleus, M/j are the nuclear Gamow- Teller 
and Fermi type transition probabilities [ p^ . 

In a nucleus as heavy as Pb the distortion of the out- 
going lepton wavefunction due to the Coulomb field of 
the daughter nucleus becomes large and affects the in- 
tegrated cross section considerably. In our treatment of 
this effect we follow the findings of ref. In ref. Q 
it is found that the "Effective Momentum Approxima- 
tion" (EMA) works well for high energy neutrinos. This 
approximation consists in using an effective momentum 

pf^ = ^fE^~^ where Eefj ^ E - VciO) (Fc(0) 

is the Coulomb potential at the origin) in calculating 
the angle integrated cross section and multiplying eq.(||) 
by {pf^/pi)^- It is also shown that the Modified EMA 
(MEMA) works better than EMA for i^^ of low and high 
energies. In this approximation eq.(|l|) is multiplied by 
pI^ ^ Eef f /piEi. We use therefore this method in all our 
calculations of the {v^,^~) cross sections. In the case 
of the (i'f,T&~) process, the situation is somewhat more 
complicated. The Fermi function works only for very low 



energies, namely Ef. < 10 MeV (where PeR ^ 1, i? is the 
nuclear radius), whereas the EMA seems to be a good ap- 
proximation for most energies of the outgoing electrons 
]l^. As in ref. |l|, for i/g [|o|, we treat Coulomb correc- 
tions by interpolating between the Fermi function at low 
electron energies and the EMA approximation at high 
lepton energies. 

To get flux-averaged cross sections it is necessary to 
convolute (|^) by the neutrino flux f{E„), that is 

{a)f = / dE,a{E,)f{E,), (2) 
Je„ 

Eq being the threshold energy. The choice of /{E^) de- 
pends on the neutrino source and can be taken for ex- 
ample equal to the supernova neutrino energy spectrum 
given by transport codes or the neutrino fluxes produced 
by a beam dump. 

The nuclear structure model used to evaluate the tran- 
sition probabilities Mjs in (|l|) is the charge- exchange 
Random-Phase- Approximation (RPA) . The details of the 
approach can be found in ||2l[. The calculations we 
present have been obtained in a self-consistent approach: 
the HP single-particles energies and wavefunctions as well 
as the residual particle-hole interaction are derived from 
the same effective forces, namely the SIII ||2^ and SGII 
psf Skyrme forces. We have found that the model con- 
figuration space used is large enough for the Ikeda and 
Fermi sum rule to be satisfied as well as the non energy- 
weigthed and energy-weighted sum rules for the forbid- 
den transitions ||2^. The GT strength distribution we 
have obtained is peaked at 19.2 MeV, in agreement with 
the experimental value. This main peak exhausts about 
60% of the Ikeda sum rule. The IAS results at 18.4 MeV 
and this value compares again well with the experimen- 
tal finding (18.8 MeV). Apart from these two resonances 
and the spin-dipole, the experimental knowledge about 
states of higher multipolarity is rather poor. The recent 
experiment of Ref. P5[ shows that isovector monopole 
strength exists in '^^^Bi between 30 and 45 MeV and in 
the present calculation wc find some strength in the same 
energy region. 

In figs.l and 2 we show the non fiux-averaged 
2"8p6(j/e,e-)208Bi and Pb{iy^, fi-)^°^ Bi inclusive 
cross sections as a function of the neutrino energy, for 
a mesh of energies, namely AE = 2.5 MeV for E,y^ and 
AE = 5.0 MeV for i?^^. The dashed fine in fig.l shows 
the cross section obtained when only the Fermi function 
is used to include the Coulomb corrections. The results 
shown have been obtained with the SIII force, but we 
have found that with the SGII force we get quite sim- 
ilar results. All the multipolarities with J < 6 are in- 
cluded. We have checked that the contribution coming 
from J = 7 is small. (Note that, for higher multipolari- 
ties, a mean field description, neglecting the particle-hole 
residual interaction, can be used to evaluate the transi- 
tion probabilities (fl)). In the calculations we present the 
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axial vector coupling constant has been taken equal to 
1.26. Note that the use of an effective ga to take into 
account the problem of the "missing" GT strength will 
reduce the reaction cross section by 10 — 15% as it was 
already discussed in ]l6| |. 

Figure 3 shows the contribution of the different multi- 
polarities to the total cross section (fig.l), for the imping- 
ing neutrino energies E^^ = 15,30,50 MeV, which are 
characteristic average energies for supernova neutrinos. 
When Ei,^ = 15 MeV (fig. 3, up), cfi,^ is dominated by 
the allowed Gamow- Teller {J^ = 1+) transition. As the 
neutrino energy increases (fig. 3, middle), the allowed IAS 
and other forbidden transitions start to contribute sig- 
nificantly. Finally, when E^^ = 50 MeV (fig. 3, bottom), 
the GT and IAS transitions are not dominating at all, 
the cross section is being spread over many multipolari- 
ties. These results suggest that r-process nucleosynthesis 
calculations such as 1^ , which include neutrino- nuclei re- 
actions, should take into account forbidden transitions. 



This may be even more important if Vj- , v^^ 
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tions occur, because in this case electron neutrino may 
have a higher average energy than it is usually expected 
from current supernovae models. 

Let us now come to the comparison with other avail- 
able calculations. Table 1 shows our flux-averaged cross 
sections, in comparison with those of refs. |l^Jl^. The 
low-energy neutrino flux is given by a Fermi-Dirac spec- 
trum [0,[5) 

1 El 
■^^^"^ " c(a)r3 exp [(E^/T) ~ a] + 1 

where T, a are fitted to numerical spectra and c{a) nor- 
malizes the spectrum to unit flux. The values of the 
parameters T and a have been chosen to be able to com- 
pare our results with those of ||l^,|l5| . As we can see from 
table 1, our predictions are in close agreement (the dif- 
ference is at most 20-30%) with The results of ||l|] 
have been obtained in a CRPA approach. A variation of 
20 — 30% is actually to be expected for calculations based 
on the same approach but using different parametriza- 
tion (for example for single particle wavefunctions and 
effective particle- hole interaction), because of the sensi- 
tivity of the flux-averaged cross sections to the detailed 
strength distributions |l^ , as we will discuss further. On 
the contrary, our results and those of ref. [p^ present sig- 
nificant differences with those of jl^ , obtained using the 
allowed approximation and including the IAS, the GT 
and the first-forbidden contributions treated on the basis 
of the Goldhaber- Teller model. 

We have checked that the differences do not come from 
the higher order forbidden transitions which are not in- 
cluded in the calculations of |l^]. The three calculations 
satisfy the same constraints, namely they reproduce the 
centroid of the resonances and satisfy the sum-rules. 

We believe that the significant differences (by a factor 
of 2) with 0| may have two origins. The first possi- 



ble origin might be the way the Coulomb corrections are 
treated. In |l^, the Coulomb distortion of the outgoing 
electron wave function was taken into account by multi- 
plying the cross section (|l|) by a Fermi function. In order 
to see the effect of using only the Fermi function instead 
of making an interpolation between the Fermi function 
and the EMA approximation, we have calculated the re- 
action cross sections using these two possible corrections. 
As figure 4 shows, the two cross sections have a quite 
different behaviour as a function of the neutrino energy 
so that this difference on the flux-averaged cross section 
may vary according to the particular neutrino flux con- 
sidered. To get a quantitative idea of the variation, we 
have calculated the flux-averaged cross sections by con- 
voluting the two curves of fig. 4 with (|^). If we use the 
Fermi function only, the reaction cross sections increase, 
on average, by 50%. 

The second possible origin of the discrepancies between 
our work, ||l^ and |l^ might be the sensitivity of the 
flux-averaged cross sections to the detailed strength dis- 
tributions in ^'^^Bi. In fact, it has already been dis- 
cussed in ]T^ , that for low-energy neutrinos, the flux- 
averaged cross sections are very sensitive to the energy 
of the excited states in the final nucleus. The reason is 
twofold. First, due to the small electron mass, the non 
flux-averaged cross section (Q) scales as the square of en- 
ergy of the states. Second, the energy dependence of the 
neutrino flux may emphasize differences in the non flux- 
averaged cross sections due to variations in the energy of 
the states. As it was discussed in [|l^, these two effects 
may modify the flux-averaged reaction cross sections by 
20 - 30%. 

To complete our comparison with the calculations of 
p^ , we have calculated two more flux-averaged cross 
sections, using the neutrino fluxes of both i/^ coming 
from the DIF of 7r+ and ly^ coming from the DAR of 
The neutrino fluxes f{E^) were taken from [p6| . 
These neutrino fluxes have been used in the recent ex- 



periments 

performed by the LSND and KARMEN collaborations. 
The DAR{i'e,e~) cross section calculated is aoAR — 
44.39 • 10"-*'' cm^ which is very close to 36.2 • 10-'^° err? 
obtained in [Tsl 
(7DIF = 399.2 



On the contrary, our DIF{v^j^ 
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cm 



^: whereas the one of ItK 



115T0~ cm . We believe that some of the disagreement 
may come from differences in the strength distributions 
of the high order (higher than 2) forbidden transitions. 
In fact, contrary to the reactions of neutrinos on light 
nuclei such as carbon, where these states contribute only 
by 20% to the total DIF cross section, their contribution 
represents 65% of the total cross section when the nucleus 
is as heavy as lead. 

Let us finally discuss the two possible neutrino oscil- 
lation u^^Vr — » Ve signals based on the spectroscopic 
properties of ^^^Bi excited in the CC reaction that have 
been proposed recently. In ||l^ , it was shown that the 2- 
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neutron events associated with the deexcitation of "^^^Bi 
are very sensitive to the mean electron neutrino energy. 
This signal relies on the fact that most of the IAS, GT 
and first-forbidden strength distributions are above the 
2n emission threshold (14.98 MeV) in '^'^^Bi. Our re- 
sults show that not only the allowed and spin-dipole 
strengths are above this threshold, but also a fraction 
of the strength distributions associated with other for- 
bidden transitions (fig. 3) will contribute to the 2n decay. 
All the arguments given in are based on the statis- 
tical calculations of \n and 2n decays. The direct In 
emission represents about 50% of the total width in the 
case of the IAS, and 5 - 10% in the case of the GT [|ll. 

In ||l5| , it was pointed out that the energy distribution 
of the neutrons in the In events should form a peak at 
low energy, more or less pronounced according to the oc- 
curence or absence of oscillations. This peak comes from 
the GT strength distribution at around 7.6 MeV which is 
located above the \n threshold emission at 6.9 MeV . Our 
GT distribution also shows a peak at around 7.5 MeV. 
We have checked that its location is not sensitive to the 
choice of the effective forces used. Still one should be 
careful about conclusions, because predictions of differ- 
ent models about the energy location and strength of that 
peak are at variance. 

In summary, we have presented the non flux-averaged 
'^°^Ph{ve,e-'f°^Bi and ^osp^^^^^ ^-^208^^ reaction 
cross sections, calculated in a self-consistent charge- 
exchange Random-Phase- Approximation with Skyrme 
effective forces. These predictions can be employed for 
very different purposes, such as for the interpretation 
of the recent experiments on neutrino oscillations per- 
formed by the LSND collaboration (where reactions in- 
duced by neutrinos on lead contribute significantly to 
the background) and to evaluate the feasibility of future 
projects in which lead should be used as detector for su- 
pernova neutrinos. We have emphasized that forbidden 
transitions contribute significantly to the neutrino-nuclei 
reaction cross sections even at the "astrophysical neu- 
trino energies" and they should be included in present r- 
process nucleosynthesis calculations. We have discussed 
the present status on the theoretical predictions on the 
reaction cross sections for the Ue having typical energies 
from present models on supernovae. If on one hand our 
calculations agree with those of ref. |l^, which are also 
based on RPA; on the other hand, they both significantly 
disagree with those of ref. ||l^. We point out that the 
origin of the discrepancy might be mainly the different 
treatment of Coulomb corrections, but also the sensitiv- 
ity of the reaction cross sections to the detailed energy 
spectrum of the final nucleus. We have also compared 
our flux-avcragcd reaction cross sections with v^^ com- 
ing from the DIF of tt^ and with coming from the 
DAR of ^JL+, with the ones of |15| ]. As expected, the DAR 
cross sections are very close. On the contrary our DIF 



have pointed out that the two predictions may differ be- 
cause of differences in the strength distributions of for- 
bidden transitions of high multipolarity which represent 
the main contribution in reactions of neutrinos on nuclei 
as heavy as lead. Finally, we have discussed our results 
in relation with recently proposed signals to measure su- 
pernova neutrino oscillations. 

This work was supported by the US-Israel Binational 
Science Foundation. 
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FIG. 1. DifTerential ^''**P6(i^e, e~)2°®Bi cross section as a 
function of electron neutrino energy for a mesh of energies 
(Ai?^^ = 2.5 MeV). As far as the treatment of the Coulomb 
distortion is concerned, the results are obtained by interpolat- 
ing between the Formi function, good at low electron energies 
and the Modified Effective Momentum Approximation, good 
at high outgoing electron energies. The dashed line shows the 
result obtained if only a Fermi function is used. 



TABLE I. Flux-averaged cross sections (10 



2) ob- 



tained by convoluting the inclusive cross sections of fig.l by 
a Fermi-Dirac spectrum (3) for neutrinos emitted in a super- 
nova explosion. Different temperatures T and a values are 
considered. The results of recent calculations are shown for 
comparison. 
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FIG. 2. Differential '^"^ Pb{Uf„ Bi cross section, ob- 

tained with the MEMA approximation, as a function of muon 
neutrino energy for a mesh of energies (AE^ = 5.0 MeV). 
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FIG. 3. Contribution of the different multipolarities to the 
differential ^°**Pfe(i/e, e~)^°®Bi cross section (10"*" cm^) of 
fig.l for E^^ = 15 MeV (up), 30 MeV (middle), 50 MeV 
(bottom). 
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